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Summary. Dynamical scaling and ageing in disordered systems far from equilib- 
rium is reviewed. Particular attention is devoted to the question to what extent a 
recently introduced generalization of dynamical scaling to local scale-invariance can 
describe data for either non-glassy systems quenched to below T c or else for spin 
glasses at criticality. The dependence of the scaling behaviour on the distribution of 
the random couplings is discussed. It is shown that finite-time corrections to scaling 
can become quite sizable in these systems. Numerically determined ageing quantities 
are confronted with available experimental results. 



1 Introduction 

Understanding cooperative phenomena far from equilibrium poses one of the most 
challenging research problems of present-day many-body physics. At the same time, 
the practical handling of many of these materials has been pushed to great so- 
phistication, and a lot of practical knowledge about them exists since prehistoric 
times. Paradigmatic examples of such system are glasses. In many cases, they are 
made by rapidly cooling ('quenching') a molten liquid to below some characteristic 
temperature-threshold. If this cooling happens rapidly enough, normal crystalliza- 
tion no longer takes place and the material remains in some non-equilibrium state. 
These non-equilibrium states may at first and even second sight look very stationary 
- everyone has probably seen in archaeological museums intact specimens of Roman 
glass or even older tools from the Paleolithic or old-stone-age - after all, obsidian or 
fire-stone is a quenched volcanic melt. But since the material is not at equilibrium, 
at least in principle it is possible (and it does happen very often in practice) that 
over time the properties of the material change - in other words, the material ages. 1 
The properties of such non-equilibrium systems depend on the time - their age - 
since they were brought out of equilibrium and this is colloquially referred to as 
ageing behaviour. 

1 Recall that physical ageing as it is understood here comes from reversible micro- 
scopic processes, whereas chemical or biological ageing may come from the action 
of essentially irreversible (bio-) chemical processes. 
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Although the first systematic studies of ageing were carried out in the by now 
classic experiments of Struik [116] in polymeric glasses, it has been realized in re- 
cent years that very similar phenomena already occur in non-disordered and non- 
frustrated systems which hence are also commonly referred to as ageing. In order 
to present the essence of the ageing phenomenon, we begin by considering a simple 
Ising model, made from spin variables crj = ±1 attached to each site i of a hyper- 
cubic lattice and an interaction as given by the classical hamiltonian with the usual 
nearest-neighbour interactions 

H = -J^2<7i<7j (1) 

where J > is the exchange integral and the sum extends over pairs of nearest 
neighbours. In d > 1 dimensions, this system has a phase-transition at some critical 
temperature T c > which separates a disordered phase at high temperatures with 
a single thermodynamically stable state from an ordered phase at low temperatures 
with two thermodynamically stable equilibrium states. It is well-established that 
such simple Ising models may be used to discuss uniaxial magnets, binary alloys or 
even liquid crystals. The motion of the spins is generated by coupling the model to 
a thermal bath of temperature T. A possible way of realizing this is through the 
so-called heat-bath dynamics which is defined by the stochastic rule 

ai(t + At) = ±1 with probability [1 ± tanh(/i;(t)/T)] /2 (2) 

where At is the time increment, the local time-dependent field is hi(t) = ~}2 y ^ a y(i){~t) 
and y(i) runs over the nearest neighbours of the site i. It is well-known that this 
rule satisfies detailed balance and hence the system evolves towards the equilibrium 
probability distribution P oq = Z^ 1 exp(— H/T), where Z is the canonical partition 
function [118, 93, 123]. The system is prepared at some initial temperature Ti n j far 
above the critical temperature T c > 0. The initial time t = is defined by coupling 
the system to the thermal bath at some low temperature T < T c and starting the 
dynamics. During the simulation, the temperature T is kept fixed and one observes 
the time-dependence of observables such as correlation functions or susceptibilities. 2 
Qualitatively, the behaviour of the system can be illustrated through the equi- 
librium free energies at the temperatures Ti n i and T, see figure 1. Before the quench, 
the system is at equilibrium with respect to the initial temperature T- ln i 3> T c and 
sits at the minimum of the free energy, as indicated by the black ball in figure la. 
If one would perturbe this system slightly, it would relax rapidly, i.e. with a finite 
relaxation time < r < oo, towards this unique equilibrium state. On the other 
hand, immediately after the quench the system does not yet have had the time to 
evolve but, with respect to the new equilibrium, its free energy is no longer minimal, 
see figure lb. Rather, two new local minima of the free energy appear, which corre- 
spond to the two equivalent ordered states of the system. Because of the competition 
between these two equivalent equilibrium states, the system as a whole cannot relax 
rapidly to one of them but rather undergoes a slow dynamics, with formally infi- 
nite relaxation times. Locally, each spin will be subject to the time-dependent field 

2 The chosen dynamics is such that the total average magnetization M(t) — 
(cr{t)} = X)j( cr <(^)) remains at its initial value M(0) = 0. Throughout, we shall 
always use this initial condition unless explicitly stated otherwise. 
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Fig. 1. Free energy F(M) as a function of the magnetization M of a simple ferro- 
magnet at (a) an initial high temperature T > T c before the quench and (b) after 
the quench to a low temperature T < T c . 

hi(t) coming from its neighbours and this field will tilt the balance between the two 
equivalent equilibrium states of figure lb in favour of one or the other. Physically, 
this means that the system will rapidly decompose into ordered domains and the 
slow long-time dynamics of this domain growth will be determined by the motion 
of the domain walls between these ordered domains. This slow (non-exponential) 
dynamics is the first defining property of ageing systems. 

Another aspect of this non-equilibrium dynamics (since in a spatially infinite 
system none of the equilibrium states will be reached in a finite time) becomes 
apparent if one considers a quantity like the two-time autocorrelation function of 
spins at site i at times t and s 



which by spatial translation-invariance is independent of the chosen site i. In fig- 
ure 2a data for C(t, s) plotted over against the time difference t — s are displayed 
for the three-dimensional Ising model. We see that the data depend on both t — s 
and s, hence time-translation invariance is broken and the system ages. Further, 
with increasing values of the waiting time s, the system becomes 'stiffer' and a 
plateau close to the equilibrium value C cq = M<? q develops when t — s is not too 
large before the correlations fall off rapidly when t — s — > co. Together with the 
slow dynamics mentioned above, this breaking of time-translation invariance is the 
second defining property of ageing systems. In principle, this could mean that the 
details of the dynamics of ageing systems might depend on the entire prehistory 
of the sample under study, which would make any attempt to formulate a general 
theory for such systems hopeless. However, a great simplification, due to dynamical 
scaling, is apparent in figure 2b where the same data for C(t, s), when plotted over 



C(t,s) = {ai(t)ai(s)} 



(3) 




against t/s, neatly collapse onto a single curve, if only the time s is large enough 
and also t — s is not too small. Since in domain coarsening of simple magnets one 
expects that the linear size of the ordered domains is L = L(t) ~ t 1//z when t is large 
enough and z is the dynamical exponent [20], the collapse in figure 2b means that 
C(t,s) = f(L(t)/L(s)), or in other words L(t) is the only relevant length-scale at 
time t. Dynamic scaling (although more general forms for L(t) than simple power 
laws are of course possible) is the third essential property of ageing systems. 

These three basic properties of ageing systems are also found in glassy systems. 
An important property of glasses is the presence of frustrations which prevent the 
relaxation of all local degrees of freedom. In consequence, the free-energy landscape 
of glasses can be very complex, with many local minima. The classic example for 
ageing behaviour was observed by Struik [116] in studying the mechanical properties 
of polymeric glasses which after a quench from the molten phase to low tempera- 
tures (i) relax very slowly (typical time-scale of years), (ii) show clear evidence of the 
breaking of time-translation invariance and furthermore, (iii) the experimental data 
for the time-dependent creep curves of the mechanical response can all be mapped 
onto a single master curve. Remarkably, that master curve turned out to be the same 
for materials as different as polymers such as PVC or PMMA, sugar or even metals 
like lead ! Evidently, there are universal scaling functions in ageing which exactly 
because of their universality one may hope to be able to understand theoretically. 
Returning to simple magnets, that universality in the kinetics of coarsening (with a 
non-conserved order-parameter) after a quench to T < T c is captured through the 
celebrated Allen-Cahn equation [3] which states that the velocity v of the domain 
walls which separate the ordered domains is related to the curvature K in d dimen- 
sions via the purely geometric relationship v — (d— 1)K, quite independently of any 
details of the interactions of the spins. 

The microscopic evolution which underlies the above statements about the prop- 
erties of macroscopic observables is illustrated in figure 3. In the first column, snap- 
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Fig. 3. Snapshots illustrating the growth of clusters in three variants of the kinetic 
Ising model, at times t = 25, 100 and 225 after the quench, from top to bottom. 
The left column shows the two-dimensional Glauber-Ising model without disorder 
quenched to T = 1.5 < T c and the middle column shows the random bond Ising 
model in two dimensions with e — 2 and T = 0.7 < T c . The right column shows a 
section of the three-dimensional Ising spin glass with binary disorder at T = 0.8 < 



shots of configurations of the two-dimensional Glauber-Ising model without disorder 
are shown for three different times after the quench from a fully disordered initial 
configuration into the ordered phase. Clearly, ordered domains form very rapidly and 
continue to grow up to the system size. In order to appreciate the effects of disorder, 
we next consider a random-bond Ising model defined by the classical hamiltonian 

Hdis = — ^ JijCTiCTj (4) 

where the Jjj are random variables uniformly distributed over the interval [1 — 
e/2, 1 + e/2] with < e < 2. Snapshots from this two-dimensional model are shown 
for e — 2 and T — 0.7 < T c in the middle column of figure 3. Again, detectable 
ordered domains from rapidly but with respect to the non-disordered Glauber-Ising 
model, although the clusters look broadly similar, one sees that the growth of the 
clusters proceeds much more slowly. Finally, we consider an Ising spin glass, which 
can again be defined by the hamiltonian (4) but now the random couplings Jij = ±1 
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can take both positive and negative values with equal probabilites. The right column 
in figure 3 illustrates the evolution of the three-dimensional Ising spin glass with 
binary disorder (Jjj = ±1 with equal probability) at the temperature T = 0.8 < T c . 
It is clear that the presence of both ferromagnetic and antiferromagnetic interactions 
no longer permits the formation of ordered clusters and with respect to the two 
other cases, the evolution appears considerably slower. We shall formalize this later 
on when we discuss the growth laws of the time-dependent characteristic length. 

Universality as seen in experiments [116] or manifest in the Allen-Cahn equa- 
tion may signal the presence of deeper dynamical symmetries. In this review, we 
shall describe the available evidence which points to the possibility that the well- 
established dynamical scaling found in ageing systems may be included into a larger 
group of local scale-transformations. As a necessary preparation we shall begin by 
briefly recalling the main features of dynamical scaling. The topic has been under 
intensive study, see [20, 29, 34, 49, 33, 82, 25, 46, 64, 67] for reviews. 

It is useful to discuss the properties of the coarse-grained order-parameter, de- 
noted by <f>(t,r) at time t at the location r. Most of our discussion is for the case 
that the order-parameter is non-conserved by the dynamics and that the initial state 
is totally disordered with ((j>(0, r)} — 0, unless explicitly stated otherwise. It has be- 
come common to study ageing behaviour through the two-time autocorrelation and 
(linear) autoresponse functions 



C(t,s) = (<f>(t,r)<f>(s,r)) ~s- b f c (t/s) (5) 
~ s-^Mt/s) (6) 



Sh(s, r) 

where h(s,r) is the magnetic field conjugate to (j> at time s and location r. The 
scaling behaviour is expected to apply in the so-called ageing regime where 

i, S ^> ^micro and t S ^> imicro (7) 

where t m icro is a microscopic time-scale. We emphasize the importance of the third 
condition for the validity of dynamical scaling. This makes the above qualitative 
observations drawn from figure 2 more precise. In particular, we see in figure 2 that 
for finite values of s scaling is not perfect: indeed corrections to dynamical scaling 
are detectable when f/s« 1. 

The distance of such systems from a global equilibrium state can be measured 
through the fluctuation-dissipation ratio, defined as [35] 



X(t,s):=TR(t,s)(^^j . (8) 

At equilibrium, X(t,s) — 1 from the fluctuation-dissipation theorem. One often 
considers the limit fluctuation-dissipation ratio 3 



X x := lim lim X(t,s)) = lim lim X(t,s)\ . (9) 

s^oo \t^oo / y^oo \t,s^eo \y — t/sj 

For quenches to below T c , one usually has Xoc — but for critical quenches onto 
T = T c , it has been proposed by Godreche and Luck that X x should be a uni- 
versal number [48], since it can be written as a ratio of two scaling amplitudes. 



The order of the limits is crucial, since lim^oc (lim^oc X(t, s)) = 1. 
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The universality of X^, has been thoroughly confirmed for systems relaxing towards 
equilibrium steady-states, and can be theoretically validated in a field-theoretical 
setting, see [33, 25] for recent reviews. On the other hand, we mention here that quite 
similar scaling behaviour is also found for systems containing irreversible chemical 
reactions and hence relaxing towards non-equilibrium steady-states. In that case, 
the definitions of X(t, s) and of Xoo have to be reconsidered, see [67] for a recent 
review. 

Furthermore, in writing eqs. (5,6) it was tacitly assumed that the scaling derives 
from the algebraic time-dependence of a single characteristic length-scale growing 
algebraically with time Lit) ~ t 1 ^ which measures the linear size of correlated or 
ordered clusters and where z is the dynamic exponent. That is indeed the case for 
simple magnets whose kinetics may be described in terms of a simple Ising model, 
e.g. with Glauber-type dynamics, where it can even be shown that z = 2 [20, 21]. 
Then the above forms (5,6) define the non-equilibrium exponents a and b and the 
scaling functions fc(y) and fn(y). For large arguments y — > co, one generically 
expects 

fc(y) ~ y' Xc/z , f R (y)~y- XR/z (10) 

where Ac and \r, respectively, are known as autocorrelation [44, 73] and autore- 
sponse exponents [102]. 

While in non-disordered magnets with short-ranged initial conditions one usually 
has Ac = Ajj, this is not necessarily so if either of these conditions is relaxed. From a 
field-theoretical point of view it is known that for a non-conserved order-parameter 
the calculation of Xc,r requires an independent renormalization and hence one can- 
not expect to find a scaling relation between these and equilibrium exponents (in- 
cluding z) [76]. On the other hand, the values of the exponents a and b are known. 
For quenches to T = T c , the relevant length-scale is set by the time-dependent cor- 
relation length L(t) ~ £(i) ~ t 1 ^ and this leads to a = b = (d — 2 + rj)/z, where 
77 is a standard equilibrium exponent. For quenches into the ordered phase T <T C , 
one usually observes simple scaling of C(t,s) = fc(t/s), hence b — 0. 4 The value 
of a depends on whether the equilibrium correlator is short- or long-ranged, respec- 

Table 1. Values of the non-equilibrium exponents a, b and z for non-conserved 
ferromagnets with T c > and a vanishing initial magnetization. The non-trivial 
critical-point value z c is model-dependent. 





a 




b 


z 


Class 


T = T C 


(d- 


2 + n)/z 


(d-2 + r])/z 


Zc 


L 


T <T C 


(d- 


2 + v )/z 





2 


L 




1/z 







2 


S 



4 This needs no longer be the case when the ageing close to a free surface is con- 
sidered [10]. 
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tively. These may be referred to as classes S and L, respectively, and one has, see 
e.g. [29, 56, 60] 

n i \ fe~ l?,|/e J class S _ J 1/z , . 

° c,lrJ ~ \ | r |-(d-2+r,) =>• | c i ass L ^ a ~ \ [d-2 + n)/z [ ' 

Examples for short-ranged models (class S) include the Ising or Potts models in 
d > 1 dimensions (and T < T c ), while all systems quenched to criticality, or the 
spherical model or the two-dimensional XY model below the Kosterlitz-Thouless 
transition are examples for long-ranged systems (class L). 

The values of these exponents are again collected in table 1. They have been 
extensively confirmed in many numerical studies, see [56, 60, 30, 1, 2, 89], and are 
reproduced in all known analytically solvable models. 

In section 2, we recall some of the main facts related to an extension of dynamical 
scaling towards a local scale-invariance (LSI), first for a dynamical exponent z = 
2 and later on possible extensions for z / 2 and review its application to non- 
disordered systems. In section 3 we discuss the scaling of responses and correlators 
in disordered, but non-glassy systems and shall present the available evidence that 
LSI should be extendable to this class of systems. In section 4 we consider the critical 
Ising spin glass, present evidence that its non-equilibrium properties depend on the 
distribution of the coupling constants and show to what extent the data can be 
explained in terms of LSI. We summarize in section 5. 



2 Local scale-invariance without disorder 

In equilibrium critical phenomena, it is well-known that the standard scale-invariance 
can, under quite weak conditions, be extended to a conformal invariance. Roughly, a 
conformal transformation is a scale-transformation r i— » br with a space-dependent 
rescaling factor b = b(r) (such that angles are kept unchanged). In particular, in two 
dimensions conformal invariance allows to derive from the representation theory of 
the conformal (Virasoro) algebra the possible values of the critical exponents, to set 
up a list of possible universality classes, to calculate explicitly all n-point correlation 
functions and so on, see [14, 28, 45, 53]. One might wonder whether a similar exten- 
sion might be possible at least in some instances of dynamical scaling and further 
ask whether response functions or correlation functions might be found from their 
covariance under some generalized dynamical scaling with a space-time- dependent 
reseating factor b — b(t,r) [52, 55] ? 

We shall first consider the case of phase-ordering where L(t) ~ t 1 ^ with z = 2 is 
known [20] and later describe how this might be generalized to more general values 
of z 7^ 2. In particular we enquire what can be said about the scaling functions 
fc,n{y) in a model-independent way. 

A useful starting point is to consider the symmetries of the free diffusion (or free 
Schrodinger) equation 

2Md t <t> = A<f> (12) 

where A — V • V is the spatial laplacian and the 'mass' M plays the role of a 
kinetic coefficient. Indeed, it was already shown by Lie more than a century ago 
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that this equation has more symmetries than the obvious translation- and rotation- 
invariances. Consider the so-called Schrodinger-group defined through the space-time 
transformations 

tt-^t'= at ~^~ @ ■ r |_> r ' — ^ r + vt _t a a g _ q _ i MgN 
ft + 5 7t + <5 

where a, /3, 7, 5, v, a are real (vector) parameters and R is a rotation matrix in 
d spatial dimensions. The group acts projectively on a solution <f> of the diffusion 
equation through (t, r) 1— » g(t, r), <f> 1— » T 9 

(T S «A) (t, r) = /^(t, r)) ^(t, r)) (14) 

where p is an element of the Schrodinger group and the companion function reads 
[95, 100] 



/ 9 (t,r) = ( 7 t + ,Sr d/2 e X p 



M -yr 2 + 2 R r ■ (7a - 5v) + 7a 2 - tSv 2 + 2 7 a • v 



2 yt + 5 

(15) 

It is then natural to include also arbitrary phase-shifts of the wave function <f> within 
the Schrodinger group Sch(d). In what follows, we denote by sch d the Lie algebra 
of Sch(d). The Schrodinger group so defined is the largest group which maps any 
solution of the free Schrodinger equation (with M fixed) onto another solution. For 
a simple illustration, consider the case d = 1 and define the Schrodinger operator 

S := 2M X_i - Y 2 1/2 . (16) 

The Schrodinger Lie algebra scfy 1 = (X-1,0,1 , Y_ 1 i,Mo) is spanned by the in- 
finitesimal generators of temporal and spatial translations (X-i, Y_\ / 2 ), Galilei- 
transformations (Vi/2), phase shifts (M ), space-time dilatations with z — 2 (X ) 
and so-called special transformations (Xi). Explicitly, the generators read [52] 

v ,»+i a n + l n „ n(n + l) n _x 2 x, „ 

X n = -t ^ d t — * rd r -Alt r - -(ra + 1)£ 

2 4 2 

y m = -t m+1/2 d r - (m + ^ Mt m - 1/2 r (17) 
M n = -Mi™ 

Here a; is the scaling dimension and M is the mass of the scaling operator <j> on 
which these generators act. The non-vanishing commutation relations are 

[X„,X n >] = (n - n')X n+n i , [X n ,Y m ] = (J^-rn)Y n+m 

[X n ,M n ,] = -n'M n+n , , [y m ,y m ,] = (m-m')M m+m , (18) 

The invariance of the diffusion equation under the action of sct) 1 is now seen from 
the following commutators which follow from the explicit form (17) 

[<S,X_i] = [5,y ±1/2 ] = [5, M ] = 
[S,X ] = -S , [<S,Xi] = -2tS - (2x - 1)M (19) 

Therefore, for any solution <j> of the Schrodinger equation S<p = with scaling di- 
mension x = 1/2, the infinitesimally transformed solution X(f> with X € sci) 1 also 
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satisfies the Schrddinger equation SX<j> = [80, 95, 51]. For applications to ageing, 
we must consider to so-called ageing algebra age 1 = (Xo,i, Y_i i, Mo) C st\) 1 (with- 
out time-translations) which is a true subalgebra of sd) 1 . Extensions to d > 1 are 
straightforward. 

What is the usefulness of knowing dynamical symmetries of free, simple diffusion 
for the understanding of non-equilibrium kinetics ? One way of setting up the prob- 
lem would be to write down a stochastic Langevin equation for the order-parameter. 
The simplest case is usually considered to be a dynamics without macroscopic con- 
servation laws (model A), where one would have [71] 

™%-»-*& + , (») 

where V is the Ginzburg-Landau potential and ij is a gaussian noise which describes 
the coupling to an external heat-bath and the initial distribution of <f>. At first sight, 
there appear to be no non-trivial symmetries, since the noise term r\ is incompatible 
with a Galilei-invariant equation (20). To understand this physically, consider a 
magnet which is at rest with respect to a homogeneous heat-bath at temperature T. 
If the magnet is moved with a constant velocity with respect to the heat-bath, the 
effective temperature will now appear to be direction-dependent, and the heat-bath 
is no longer homogeneous. However, that difficulty can be avoided as follows [103]: 
split the Langevin equation into a 'deterministic' part with non-trivial symmetries 
and a 'noise' part and then show using these symmetries that all averages can be 
reduced exactly to averages within the deterministic, noiseless theory. Technically, 
one first constructs in the standard fashion (Janssen-de Dominicis procedure) [36, 77] 
the associated stochastic field-theory with action J[cf>, <j>] where <f> is the response field 
associated to the order-parameter (j>. Second, decompose the action into two parts 

j[«M] = j [«M + j 6 [£] (21) 

where 

J o [<£,0]=/ dtdr 4>( 2Md t <t> - A4> + ^- ) (22) 



contains the terms coming from the 'deterministic' part of the Langevin equation 
(V is the self-interacting 'potential') whereas 

J b [4>] = -T f AtAr4>{t,rf-\[ Ardr' 0(O,r)o(r - r')4(0,r') (23) 

J«. + xM d ^ JR 2d 

contains the terms coming either from the noisy coupling to the heat bath or else 
from the average over the disordered initial conditions [77]. It was assumed here 
that (0(0, r)) = and a(r) denotes the initial two-point correlator 

a(r) := C(0, 0; r + r', r') = (0(0, r + r')4>(0, r')> = a(-r) (24) 

while the last relation follows from spatial translation-invariance which we shall 
admit throughout. Averages are of course calculated from the functional integral 

(-4) = |w^,^]e- JWl . (25) 

It is instructive to consider briefly the case of a free field, where V = 0. Variation 
of (21) with respect to 4> and <f>, respectively, then leads to the equations of motion 
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2Md t (f> = A<t> + T4> , -2Md t </> = A<j>. (26) 

The first one of those might be viewed as a Langevin equation if is interpreted 
as a noise. Comparison of the two equations of motion (26) shows that if the order- 
parameter <f> is characterized by the 'mass' M (which by physical convention is 
positive) , then the associated response field <j) is characterized by the negative mass 
— M. This characterization remains valid beyond free fields. 

We now concentrate on actions Jo [<j>, 4>] which are Galilei-invariant. This means 
that if (.)o denotes the averages calculated only with the action Jo, the Bargman 
superselection rules [7] 

~ 8 n , m (27) 



hold true. It follows that both response and correlation functions can be exactly 
expressed in terms of averages with respect to the deterministic part alone. For 
example (we suppress for notational simplicity the spatial coordinates) [103] 



Sh(s) 



{4>(W(s)) = (WMs) e~ Jbl * ] ) o = {<KW*)) (28) 



where the 'noise' part of the action was included in the observable and the Bargman 
superselection rule (27) was used. In other words, the two-time response function 
R(t,s) = Ro(t,s) does not depend explicitly on the 'noise' at all ! The correlation 
function is reduced similarly [103] 

C(t, s;r)=T f dudR U{t, r + y)<j>{s, y)4>{u, R) 2 ) 
+ \ [ dRdR' a(R-R') U(t, r + y)<j>{s,y)${Q,R)4>{Q,R!)) . (29) 

1 J R 2d \ I 

Only terms which depend explicitly on the 'noise' remain - recall the vanishing of the 
'noiseless' two-point function {<f>(t)4>(s))o = because of the Bargman superselection 
rule. Using mathematical methods explained in [57] and restricting to a disordered 
initial state a(r) — aoS(r), the resulting three-point function can be calculated. It is 
satisfying that Ac = A_r follows [103], in agreement with an earlier derivation [20]. 

Therefore, the dynamical symmetries of non-equilibrium kinetics are character- 
ized by the 'deterministic' part of Langevin equation. This result only ^iepends on 
the Galilei-invariance of the 'deterministic' part (provided that (j> and <j) transform 
projectively through eq. (14)). Clearly, linear equations are Galilei-invariant in this 
sense. On the other hand, the question is considerably more complicated for non- 
linear equations. For the purposes of this review it is sufficient to state that deter- 
ministic non-linear diffusion/Schrodinger equations with age 1 or sci) 1 as a dynamical 
symmetry and sufficiently general to be applicable to phase-ordering kinetics have 
been explicitly constructed [115, 9]. In particular, it has been shown that indepen- 
dently of the form of the potential the same representation of age d or scf) d applies, 
in agreement with universality. We refer to the literature for the details. 

Since all quantities of interest will reduce to some kind of response function, 
one may calculate them from the requirement that they transform covariantly under 
the action ageing subgroup (with Lie algebra age d ) obtained from the Schrodinger 
group when leaving out time-translations. We shall concentrate here on the two-time 
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autoresponse function R(t,s) = (4>(t)4>(s)} built from so-called quasiprimary [14] 
scaling operators <j> and 4> which transform according to the generators of the ageing 
algebra. Then the requirement of covariance of R reduces to the two conditions 
XoR(t,s) = XiR(t, s) — 0. Since time-translations are not included in the ageing 
group, the generators X n can be generalized from (17) to the following form 

X n = -r+^t-^^yVrdr- (nH ^ 1)n A4r-V-|(n+l)r-£nr ; n>0 (30) 

where £ is a new quantum number associated with the field <f> on which the generators 
X n act. Consequently, a quasiprimary scaling operator which transforms covariantly 
under the ageing algebra is characterized by the triplet (x,£,M). The last term in 
(30) can only be present for systems out of an equilibrium state. In particular, the 
requirement of time-translation invariance and [Xi,X_i] = 2Xo lead to £ = and 
any representation of scf) d must have £ = 0. In addition, the last commutator in 
(19) must be replaced by [<S, X\] — —2tS — (2x + 2£ — 1)Mq. The meaning of £ can 
be understood by integrating the generators X n in order to recover the finite, non- 
infinitesimal transformations. Then it follows that if £ 7^ 0, lattice observables are 
related to the quasiprimary fields $ of ageing invariance according to lattice (i) >- » 
a~ x P<P(t) with an additional time-dependent factor [65] and where a is a lattice 
constant. Rather than being an exotic exception, this extension of the Schrodinger 
transformations seems to occur quite generically, the best-known examples being the 
one-dimensional Glauber-Ising model at zero temperature. For further examples, see 
table 3 below. Finally, solving the two differential equations for R gives the explicit 
form of R(t,s), see (34) below. 

While this discussion was carried out explicitly for the case 2 = 2, it is tempting 
to try and generalize this idea to more general values of 2. In this way, the notion 
of local scale-transformation has been introduced, which is based on the following 
main assumptions [55]. 

1. In principle, the following conformal time-transformations should be included 

"- > ' ,= £ 7Z7 i "*-/*Y=l (31) 
■yt + d 

For applications to ageing, however, time-translations generated by /3 must be 
left out (generalizing the restriction sct) d — > agc d ). 

2. The generator Xo of scale-transformations is 

X = -tdt - -rd r - - (32) 

2 2 

where x is the scaling dimension of the quasi-primary operator on which Xo is 
supposed to act. Physically, this implies that there is a single relevant length 
scale L(t) ~ t 1/z . 

3. Spatial translation-invariance is required. 

Generators for infinitesimal local scale-transformations have been explicitly con- 
structed [55]. In almost all cases, however, it turned out that terms containing cer- 
tain fractional derivatives are needed. Skipping over the details of their construction 
which can be found in the literature, it is sufficient here to state that one can show 
that the linear equation S z <t> = with the generalized Schrodinger operator 
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S z = -\dt+(VlY 12 (33) 

where A is a constant, satisfies with the generators of local scale-transformations 
commutation relations quite similar to eq. (19) which establishes local scale- invariancc 
as a dynamical symmetry of that linear equation [55, 110, 12]. Attempts to sup- 
press the fractional terms in the generators of local scale-transformations lead ei- 
ther to A — > or to A — > oo (there is an exception if z = 1 [55]). Explicit tests 
have recently been performed for z — 4 in models of kinetic growth described by 
the Mullins- Herring equation [110] and the kinetic spherical model with conserved 
order-parameter [12]. Further tests in long-ranged models where z can be continu- 
ously tuned through a control parameter are currently being carried out [13]. On the 
other hand, a systematic derivation of an analogue of the Bargman superselection 
rules for z 7^ 2 is still lacking and indeed constitutes one of the most important open 
questions to be overcome for a further development of the theory. 

Local scale-invariance (LSI) assumes in particular that the two-time response 
functions transform covariantly under these local scale-transformations, hence XoR = 
XiR = 0. This leads to the following prediction for the autoresponse [55, 103, 62, 65] 

R(t,s)^{^(s)) = s- 1 - a f R (t,s) , fH(y)=foy 1+a '- Xa/z (y-ir 1 - a ' (34) 

where the exponents a,a,\n/z are related to x,£,x, £ and fo is a normalization 
constant. 5 Spatio-temporal responses can be found similarly, with the result 

= R(t,s)4>(\r\(t-s)- 1/z ) (35) 

fi=0 v ' 

where <P(u) — exp(— ^Mu 2 ) if z = 2. For z / 2, the function <P(u) has to be found 
from a known fractional differential equation [55, 110, 12]. 

Starting with [54], the prediction (34) has been reproduced in many different spin 
systems and we list examples quenched to below criticality in table 2 and quenched 
to the critical point in table 3. For T < T c , it is found empirically that a — a' in all 
examples considered so far. 6 We point out that agreement with local scale-invariance 
eq. (34) is not only obtained for systems where the dynamical exponent is z — 2, 
but that rather there exist quite a few examples where z can become considerably 
larger or smaller than 2. It must be remembered, however, that the above derivation 
of (34) for a stochastic Langevin equation has for the time being only been carried 
out for z = 2 and the justification of XoR = X\R = remains an open problem for 
z 7^ 2 in general although the result eq. (34) seems to work remarkably well. It is 
non-trivial that a relatively simple extension of dynamical scaling should be capable 
of making predictions which can be reproduced in physically quite different systems. 

A few comments are still needed: (i) for the XY model in the spin-wave approx- 
imation (table 2), eq. (34) holds for the response of the angular variable cj> — <j>(t, r) 
which is related to the XY spin through S = (cos (j>, sin <f>). Magnetic responses have 

5 We point out that the prediction (34) as well as the explicit form (30) of X n , 
valid for 2 = 2, assume that the mean order-parameter (4>(0, r)} = mo = at the 
initial moment when the quench to T < T c or T = T c is made. 

6 Recall that Mtrm(S, s) in direct space is not very sensitive to a — a' and that most 
numerical studies were carried out in this setting. A considerably more sensitive 
test looks at Mtrm in momentum space [104]. 



R{t,s;r) = 



6W,r)) 



Sh(s,0) 
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Table 2. Magnetic systems quenched into the coexistence phase (T < T c ) which 
satisfy (34) with the exponents a = a' and \r. d is the spatial dimension and 
the numbers in brackets estimate the numerical uncertainty in the last digit (s). 
In the spherical model, long-range initial conditions are included and in the long- 
range spherical model the exchange couplings decay as J r ~ \r\~ d ~ a . In the bond- 
disordered Ising model, the couplings are taken homogeneously from the interval 
[1 - e/2, 1 + e/2]. Then z = z(T, e) = 2 + e/T [97] and one observes roughly 1.3 < 
\ R (T,e)<1.7. 



model 


d z 


a = a' \r 




Ref. 


Ising 


2 2 

2 2 

3 2 


1/2 1.26(1) 
~ 0.5 1.24(2) 
1/2 1.60(2) 




[58] 

[89, 75] 
[58] 


Potts-3 


2 2 


0.49 1.19(3) 




[89, 75] 


Potts-8 


2 2 


0.51 1.25(1) 




[89, 75] 


XY 

XY spin wave 


3 2 
> 2 2 


0.5 1.7 
d/2-1 d 


angular response 


[2] 

[103] 


spherical 


> 2 2 


d/2 - 1 (d-a)/2 


C ini (r) ~ |r|— ° 


[94, 102] 


long-range 
spherical 


> 2 a 
< 2 a 


d/cr-1 d/2 
d/cr-1 d/2 


< CT < 2 

< a < d 


[27] 


random-bond Ising 


2 2 + e/T 


l/z(T,e) \ R (T,e) 


disordered 


[66] 



a different scaling form [17, 1]. (ii) In the critical disordered Ising model (table 3) 
one finds a logarithmic scaling form R(t, s) — (ro + ri \n(t — s))fn(t/s) [24, 111, 112] 
such that the computed /jj(j/) is consistent with (34) to one-loop order, or up to 
terms of order O(e). (iii) Finally, a two- loop calculation of the critical non-conserved 
0(n)-model does produce in 4 — e dimensions an expression for fn(y) which is in- 
compatible with (34) [23] and a similar result is anticipated in 2 + e dimensions [43] , 
although the one-loop results are still compatible [22, 23, 43]. Should one conclude 
from these studies that for T = T c the prediction (34) and by implication local 
scale-invariance can only hold approximatively ? This might well be a subtle ques- 
tion. Deviations between (34) and the field-theoretical studies typically arise when 
t/s ~ 1. However, in this region the field-theoretical results for /r(j/) do not agree 
with the ones of non-perturbative numerical studies [104] . Since the perturbative ex- 
pansion usually carried out in field-theoretical studies does not necessarily take care 
of the Galilei-invariance, it is necessary to carefully check that the truncation of the 
e-series does not introduce slight inaccuracies. Only after this has been done (for ex- 
ample by re-summing the e-series) and checked by comparing with non-perturbative 
data, meaningful quantitative statements on the scaling functions can be made, (iv) 
Throughout, it was implicitly assumed that the order-parameter vanishes initially. 
Systematic studies on what happens when this condition is relaxed are only now 
becoming available [6, 26, 43, 11]. These extensions might be particularly important 
for chemical kinetics [67] and the existing simulations for directed percolation (or the 
contact process) may well turn out to be an example where the initial non-vanishing 
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Table 3. Systems quenched to a critical point of their stationary state which satisfy 
(34) with the exponents a, a' and Xr/z. d is the spatial dimension and the numbers 
in brackets estimate the uncertainty in the last digit(s). CSM stands for the spheri- 
cal model with a conserved order-parameter, FA denotes the Frederikson- Andersen 
model, nekim is the non-equilibrium kinetic Ising model and bcp and bpcp denote 
the bosonic contact and pair-contact processes, respectively. In the bpcp dynamical 
scaling only holds along a part of the critical line. In the spherical model, long-range 
initial correlations d n i(r) ~ \r\~ d ~ a were considered. If d + a > 2, these reduce 
to short-ranged initial correlations (denoted s), but for d + a < 2 a new class L 
arises. In those models described by a Langevin equation, the simple white noise 
(ji(t, r)rj(s, r')) = 2TS(r — r')8(t — s) was used, with the only exception of the CSM. 
In the Ising spin glass, a bimodal disorder was used. 



model 


d 


a 


a' — a 


Xr/z 


Ref. 


random walk 




-1 








[35] 


OJK-model 




(d-l)/2 


-1/2 


d/4 


[16, 92, 62] 


Ising 




1 





-1/2 


1/2 






[47, 88, 59] 






2 


0.115 - 


-0.187(20) 


0.732(5) 






[104, 65] 






3 


0.506 


-0.022(5) 


1.36 






[104, 65] 


XY 


3 


0.52 





1.34(5) 


[2] 


spherical d > 2 


< 


4 


d/2 - 1 


d/2 


- a/4 - 1/2 




L 


[102] 




> 


4 


d/2 - 1 


(d 


-a)/4 + l/2 




L 


[102] 




< 


4 


d/2 - 1 





3d/4 - 1 




s 


[48] 




> 


4 


d/2 - 1 





d/2 




S 


[48] 


CSM 


> 


2 


d/4- 1 





(d + 2)/4 


[12] 


disordered Ising 


4 - 


£ 


1 _ k\fi 







0( 


s), log 


[24, 111, 112] 


FA 


> 


2 


l + d/2 


-2 


2 + d/2 






[91] 






1 


1 


-3/2 


2 






[91, 90] 


Ising spin glass 


3 


0.060(4) 


-0.76(3) 


0.38(2) 


see 


sect. 4 


[62, 63] 


contact process 




1 


-0.681 +0.270(10) 


1.76(5) 


t/s 


> 1-1 


[41, 70, 65] 




> 


4 


d/2 - 1 





d/2 + 2 






[106] 


NEKIM 


1 


-0.430(4) 





1.9(1) 


[96] 


BCP 


> 


1 


d/2 - 1 





d/2 


[8,9] 


BPCP 


> 


2 


d/2 - 1 





d/2 


a 


< ac 


[8,9] 
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value of the order-parameter influences the form of scaling functions such as /jj(y) 
[11, 41, 106, 70]. From that point of view it is surprising that eq. (34) could describe 
any part of the data of as well as it does, (v) We did not include growth 

models since their analysis from the point of view of LSI is just beginning [110, 12]. 

If z = 2, it is also possible, using eq. (29), to derive explicit predictions for the 
two-time correlation function [103, 61]. These have been tested in some exactly solv- 
able models [103, 65], the two-dimensional Ising model [61] and the two-dimensional 
g-states Potts model with q = 2, 3, 8 [89, 75]. Extensions to z = 4 have been studied 
very recently [110, 12]. 

Finally, let us mention that the prediction (35) for the space-time response has 
been verified for z — 2 through numerical simulations of two- and three-dimensional 
Ising models undergoing phase-ordering [60] and from the exact solutions of the ID 
Glauber-Ising model [47, 88, 59] at T = 0, the spherical model in d > 2 dimensions 
[48, 102] and the long-range spherical model [27, 13]. 



3 Disordered ferromagnets 

We now take up the discussion of disordered systems. Rather than going directly to 
spin glasses, it may be useful to consider first disordered, but not frustrated systems, 
such that the random exchange couplings Jjj > in models such as (4). 

We point out that a disorder in the exchange couplings, as introduced in the 
random-bond model (4), can have quite different consequences than a random site 
dilution, as described by a classical hamiltonian H = — <7X)(; j) where e, G 

{0, 1} are random variables selected according to a control parameter p. Provided 
that on average Jij — 1, the critical temperature in the presence of bond disorder 
will not change much. For a site-diluted model, however, the phase-transition should 
disappear once the non-magnetic sites (where ej = 0) will begin to percolate. Here, 
we shall mainly restrict attention to random-bond disorder and shall furthermore 
restrict to the case of quenched, that is immobile, disorder. 

The first question to be addressed concerns the existence and the detailed form 
of dynamical scaling, as expressed through the time-dependence L — L(t) of the 
linear size of the clusters. The physical picture is as follows [72], see also the middle 
column of figure 3. Initially, small clusters will form and start to grow. These early 
stages are still unaffected by the presence of the disorder which will only begin to 
be felt once the average cluster size has become of the same order as the mean 
distance between two impurities. Then the domain walls will become trapped close 
to the disorder-induced defects and their motion will slow down correspondingly. 
Rather than being described by purely geometric considerations as it can be done 
through the Allen-Cahn equation for non-disordered systems, the impurities act as 
energy barriers to the domain growth, hence the pinning centers are localized in 
energetically favourable positions, which explains the slowing-down of the kinetics. 
Phenomenologically, for a non-conserved order-parameter one still expects that the 
basic curvature-reducing mechanism is applicable, leading to [72, 84]. This yields 

dL(t) _ D(L,T) 

~dT - ~tw~ ( ' 

where the diffusion constant D — D(L, T) now depends on the domain scale L = L(t) 
and the temperature T. For a constant D, one is back to the non-disordered case, 
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with L(t) ~ t 1 / 2 . For thermally activated motion, one should have 

D(L, T) ~ D exp(-£ s /T) (37) 

where Be = Eb(L) is the barrier energy. Huse and Henly [72] have argued long ago 
that the barrier energy should depend on the domain size algebraically Eb(L) ~ 
EoL^, where the exponent tp — (2£ + d — 3)/(2 — () where £ is the domain-wall 
roughness exponent. For example, in two dimensions it is known that £ = |, hence 
ip — j. Inserting into eq. (36) leads to [72] 

T\ W /t\ _ f^r ; r«l 



which describes the qualitative change of behaviour between the initial regime and 
a second regime of slower growth for larger times. 



Table 4. Experimental results on the growth of the domain size L(t) in some disor- 
dered, but non-glassy systems in two dimensions. For Rb2Coo.60Mgo.40F4 a different 
temperature-dependent exponent ip is found close to T c . For TGS the order-parameter 
is conserved. 



material 


model 


Lit) V 


/. h\ 


Ref. 


Rb2 Coo. 60 Mgo.40 F4 


diluted Ising 
antiferromagnet 


A + B(\nt) y * 0.28(15[-K]/T) 
if T < 15[K] 


20 


[74] 


tryglycinc sulfate 
(TGS) 


ferroelectric 
conserved 


~(t-to) 1/z z~3 
~ (hit) 1 /* 1/4 


322 


[86, 87] 


Rb2Cuo.89Coo.11F4 


random-bond 
Ising 


(In*) 1 /* 0.20(5) 

T = 3.9[K] 


4.93(5) 


[114] 


ZLI 4792 (Merck) 
(liquid crystal) 


random-bond 
Ising 


(Int) 1 ^ 1/4 




[113] 



We now compare this result with the available experimental evidence collected 
in table 4. Although there is a general qualitative agreement about a a change from 
a relatively fast kinetics seen at not too late times [74, 86, 87] to slowing-down of the 
kinetics at later times, quantitatively the situation remains a little ambiguous. In 
the published work known to us, data were usually compared to a logarithmic law 
L(t) ~ (lnt) 1 /*, which rather than being an objective statement might simply come 
from the fact that the prediction (38) of Huse and Henley [72] has been around 
for quite some time. On the other hand, there is some experimental information 
which is not consistent with eq. (38), if it can be taken at face value. For example, 
in at least one system the effective exponent depends on temperature [74] while a 
constant tp = | is theoretically expected in two dimensions [72] . Indeed, it has been 
argued that in site-diluted Ising models, because of the fractal nature of the domains 
boundaries the energy barriers should rather depend logarithmically on L(t), viz. 
Eb(L) = eln(l + L) [69, 107, 97]. Inserting this into (36) leads to an algebraic law 
L(t) ~ £ 1//j where the effective 2 crosses over from z = \ for short times to 



18 



Malte Henkel and Michel Pleimling 



z = z{T,e) = 2 + e/T (39) 

for late times. In at least one experiment, data were explicitly seen to be also com- 
patible with an algebraic growth law [86, 87]. 

For a conserved order-parameter, a similar discussion can be carried out. It can 
be shown that for algebraic energy barriers the result (38) remains unchanged [72], 
whereas for logarithmic barriers one finds again algebraic growth L(t) ~ t l ^ z with 
z = z(T,e) = 3 + e/T [97, 98]. 

Theoretical work has among others studied a random-bond Ising model defined 
by the classical hamiltonian 

Hdis = - JijViVj ( 4 0) 

where the Jjj are random variables uniformly distributed over the interval [1 — 
e/2, 1 + e/2] and where a non-conserved dynamics was created via a Metropolis or 
heat-bath algorithm. For examples, estimates for L(t) extracted from the single- 
time spatio-temporal correlation function were seen to be in full agreement with 
the prediction (39) of logarithmic energy barriers [97, 98]. The same conclusions 
are obtained from the dynamical scaling behaviour of two-time quantities [66], see 
also below. It would be desirable to compare existing experimental data with the 
possibility of logarithmic energy barriers. If the discrepancy between the informa- 
tion extracted from experiments and theoretical simulations should persist, it could 
mean that the usually considered theoretical models are not realistic enough for a 
quantitative description of real disordered materials. 

We here discuss the thermoremanent magnetization Mtrm (t, s) which is ob- 
tained when quenching the system in presence of a small magnetic field. The mag- 
netic field is cut after the waiting time s and the decay of the magnetization is then 
monitored as a function of time. The thermoremanent magnetization is related to 
the response R(t, s) through 

Mtrm(M) =h [ duR{t,u), (41) 
Jo 

where ho is the amplitude of the small magnetic field. The interest in this inte- 
grated response function comes from the fact that a direct calculation of the func- 
tional derivative which defines R produces extremely noisy data, in contrast to the 
measurement of Mtrm where the noise is to a large extend smoothed out by the 
integration. 

From our experience with the scaling of Mtrm(M) m the phase-ordering of 
simple magnets, we expect the scaling behaviour [56, 60] 

Mtrm(M) =r s- a f M (t/s) +r lS - XR/z g M (t/s) , g M (y) ~ y~ XR/z . (42) 

Note that we have here included the leading correction term ~ s~ Xr > ' z guitj s) which 
can become quite sizeable in systems undergoing phase-ordering and which must 
be subtracted off before a reliable determination of the scaling function /m(?/) is 
possible. In the framework of LSI an analytical expression for /m(j/) is readily derived 
from (34) with a = a' and reads (2-F1 is a hypergeometric function) 

MV) = V' Xr/Z 2F 1 (l + a, ±* - a; ±* - a + 1; -) . (43) 
\ z z yj 
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The expected scaling behaviour (42) of the linear response has been studied 
through large-scale numerical simulations [66]. As an example, figure 4 shows data 
obtained for e = 0.5 and T = 1. This is a case where the finite-time correction to 
the leading behaviour ~ s~ a /m(£/s) is negligible, i.e. n ~ in eq. (42). Plotting 
Mtrm(us,s) over against s in a log- log plot nicely verifies the expected scaling 
behaviour and yields in addition an estimate for the ageing exponent a, see table 5. 
It is worth noting that the estimates are in complete agreement with the relation 
a(T,e) = l/z(T,e) expected for systems of class S as discussed in section 1, where 
the dynamical exponent z is given by (39). This confirms the conclusion of [97, 98] 
that the random-bond Ising model shows simple ageing and furthermore L(t) ~ t 1 ^ 
with z = z(T, e) given by eq. (39). 

We also compiled in table 5 the estimates [66] for the non-universal constants ro 
and n as well as for the exponent Xr/z which governs the asymptotic behaviour of 
the scaling function /m(j/) ~ y~ 



0.003 



0.002 - 



CO 
5^ 



DC 



0.001 - 



0.0004 




500 



Fig. 4. Thermoremanent magnetization Mtrm(us,s) of the two-dimensional 
random-bond Ising model (40) as a function of s for e = 0.5 and T — 1.0 and 
different values of y. 



The random-bond Ising model is very well suited for investigating in detail the 
possible applicability of the theory of local scale-invariance to nontrivial systems 
undergoing phase-ordering. This is of course due to the fact that the dynamical 
exponent can be continuously changed by changing the distribution of the couplings 
and/or the temperature, yielding values of z much larger than the value z — 2 
encountered in non-disordered ferromagnets undergoing phase-ordering. We have 
seen in section 2 for what models the data for Mtrm(£,s) agree with LSI, see 
table 2. Figure 5 compares for three different cases the computed scaling function 
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Table 5. Quantities describing the dynamical scaling behaviour ol the linear re- 
sponse Mtrm in the bond-disordered two-dimensional Ising model for different val- 
ues ol e and different temperatures T. The last column is obtained from equation 
(39) and should be compared with the numerical estimates for the exponent a. 



e T 


a \ R /z 


ro r 1 


1/z 


0.5 1.0 
0.8 
0.6 
0.4 


0.398(5) 0.61(1) 
0.382(4) 0.595(10) 
0.353(4) 0.58(1) 
0.310(5) 0.52(1) 


0.021(1) 
0.020(1) 
0.022(1) 
0.029(2) 0.008(1) 


0.400 
0.381 
0.353 
0.308 


1.0 1.0 
0.8 
0.6 
0.4 


0.330(5) 0.51(1) 
0.308(4) 0.49(1) 
0.277(6) 0.46(1) 
0.22(1) 0.375(10) 


0.021(1) -0.009(1) 
0.019(1) -0.006(1) 
0.020(1) -0.007(1) 
0.026(2) -0.014(3) 


0.333 
0.308 
0.273 
0.222 


2.0 1.0 
0.8 
0.6 


0.24(2) 0.33(1) 
0.22(2) 0.30(1) 
0.17(2) 0.27(1) 


0.048(2) -0.048(4) 
0.093(3) -0.042(4) 
0.194(4) -0.033(3) 


0.250 
0.222 
0.188 



/m of the thermoremanent magnetization with the LSI prediction (43). The same 
comparison was made for e = 0.5, T = 1, e = 1, T = 0.6, and e — 2, T = 1 in our 
original paper [66]. It has to be recalled that prior to this comparison the values of 
the exponents entering into (43) and also the normalization are already fixed. We 
see that the form of the scaling function (after subtraction of the above-mentioned 
finite-time correction, if needed) is perfectly described by LSI, and this is the case for 
all values of e and T investigated. This is a very remarkable result as it suggests that 
the idea of extending dynamical scaling to local, space- and time-dependent, scaling 
is capable of reproducing faithfully the linear responses of very different systems 
characterized by very different values of the dynamical exponent. 

Having looked at the evidence in favour of simple scaling of Mtrm^s), it may 
be useful here to take a slightly broader point of view. Indeed, in the discussion of 
glassy systems, people often describe the scaling of, say, the two-time correlation 
function as follows, e.g. [120] 

C(t,s) = C, t (i-a) + Cage (t,s) (44) 

such that the stationary part satisfies limt^oo C s t(t) = and furthermore 
lim t _ s ^ 00 (lim s ^ 00 C st (t — s)) = <jea, where <?ea is the Edwards- Anderson order- 
parameter for glasses. On the other hand, the ageing part is assumed to read 

Cago(M) =C (J^yJ , h(t) = h exp 

where C is a scaling function, fi is a free parameter and ho and A are constants. 
Conventionally, one refers to the case < \i < 1 as subageing, the limit case /J, — > 1 
as (full or simple) ageing and the case fi > 1 as superageing. However, it has been 
shown by Kurchan [83] that, given that the positive function C(u) decreases strictly 



A 
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ln(t/s) ln(t/s) ln(t/s) 

Fig. 5. Comparison of the numerically determined scaling functions /m(4/s) of the 
two-dimensional random-bond Ising model with the LSI prediction eq. (43) (full 
lines) for (a) e = 0.5 and T = 0.6, (b) e = 1 and T = 1 and (c) e = 2 and T = 0.6. 



monotonously with u, the case (x > 1 is incompatible with an elementary property 
of the autocorrelation function, namely that if there is a strong correlation between 
times t\ and ti > t\ and as well a strong correlation between times t% and tz > £2, 
then there must exist a strong correlation between times t\ and tz- Hence superageing 
is impossible [83]. 

There is a nice argument which explains the origin of the form (44,45) [5]. It 
relies on the observation that if C(t,s) is plotted over against t — s, a plateau is 
observed, see also figure 2a. Further, according to Zippold, Kiihn et Horner [122] 
the transition towards the scaling regime occurs at a time difference t — s ~ V" where 
< £ < 1 describes this change of behaviour. 7 Hence close to the plateau one can 
assume 

C(t, s) = f/EA + t~ a gi ((t - s)t" C ) + . . . (46) 

where a is some exponent and g\ a scaling function. Now, consider t — s = xv- and 
take the limit t — > 00 and then x ^> 1. Then, to leading order 

where b is a further exponent which can be worked out for certain model spin glasses 
[5]. Comparing eqs. (46,47), the dependences on t and on x can be separated which in 
particular leads to dln/i(i)/dt = A~ 1 t -M , where fi = f + a/b and A is a separation 
constant. Then eq. (45) follows directly [5]. In particular, this means that other 
phenomenological forms which have been discussed in the past can be eliminated. 

After these preparations we turn to a dicussion of ageing properties of the auto- 
correlation function in disordered ferromagnets. Indeed, recent studies of C(t, s) in 
bond-diluted [66] and in site-diluted [99] Ising models quenched to below the critical 
points showed surprising and at first sight unexpected features. In figure 6a we show 
data [66] for the two-dimensional random-bond Ising model plotted over against t/s. 



7 Explicitly, £ = 4/(d + 2) in the d-dimensional spherical model and £ = 4/5 in the 
spherical spin glass [122]. 
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Fig. 6. (a) Autocorrelation C(t, s) vs t/s for various waiting times s in the random- 
bond Ising model with e — 2 and T — 1. (b) Plot of C(ys, s) for the same parameters 
and y = 5, 10, 15, 20 (from top to bottom). The data can be perfectly fitted by the 
scaling form (48) with b' = 0.075, which illustrates the existence of important finite- 
time corrections, (c) The scaling function fc(t/s) obtained after subtracting off the 
finite-time corrections from the autocorrelation, again for e = 2 and T = 1, in 
agreement with simple ageing of the autocorrelation. 



In clear contrast with non-disordered simple magnets, no clear scaling is seen. At 
first sight, one might be tempted to introduce a non-vanishing ageing exponent b, 
see eq. (5), but since the curves for C(t,s) increase with s, the effective values of 
b = fecft fitted to the data will turn out to be negative [66, 99]. If true, that would 
imply that for s — » oo the autocorrelator C(t, s) would grow unboundedly, which is 
impossible. It is therefore necessary to reexamine whether the standard scaling form 
of simple ageing is applicable. Indeed, in a recent study of the randomly site diluted 
Ising model [99] the more general scaling form (44,45) was used and fitted values 
of \i in the range (i « 1.03 — 1.04 were reported, implying the so-called superageing 
behaviour [99]. Although a good collapse of the data can be obtained this way, the 
arguments raised by Kurchan [83] and quoted above make it doubtful that the con- 
clusion of [99] can be accepted at face value. Indeed, in using scaling forms such as 
(5) one has to be careful about the possibility of finite-time corrections to scaling. 
We illustrate the presence of strong corrections to scaling in figure 6b [68] for the 
case e — 2 and T = 1, for several values of y = t/s. Clearly, the data are perfectly 
described by the extended scaling form 

C(ys,s) = fc(y) ~ s~ b ' g c (y) (48) 

(hence with 6 = 0) and with b' > 0. 8 Subtracting off the leading correction term, 
yields a perfect scaling behaviour eq. (5) according to simple ageing as shown in 
figure 6c [68] , quite analogous with the standard scaling expected for simple magnets 
quenched to below their critical point. 

8 For some other values of e and of T two correction terms have to be included in 
order to describe the non-monotonic variation of C(ys, s) as a function of s, see 
[68] for details. 
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In conclusion, the random-bond Ising model apparently ages in much the same 
way as a simple magnet. The prediction of table 1 for the exponents a,b in class 
S systems hold true for this non-glassy disordered system as well, but the value 
of the dynamical exponent z = z(T, e) (and also of Ac, r) becomes dependent on 
temperature and on the distribution of the couplings. The form of the response 
function agrees with the prediction of LSI, for a large range of values of z. It would 
be very interesting to check these conclusions in other systems. 



4 Critical Ising spin glasses 

The disordered ferromagnet studied in the previous section can be viewed as in- 
termediary between the simple, non-disordered magnets and the disordered and 
highly frustrated spin glasses. In the following we discuss to what extend the re- 
sults obtained for the disordered ferromagnets can be extended to glassy systems. 
We thereby shall concentrate on the Ising spin glass, with a static Hamiltonian 
H = — Y^u,j) Ji,j a i a j- Here Oi = ±1 are the usual Ising spins and the nearest- 
neighbour couplings J; j are random variables. We shall consider three different 
distributions of the couplings: (i) the bimodal distribution with 

P B (Ji,j) = [S{Ji,j - J) + 6{Ji,j + J)]/2, (49) 

(ii) the Gaussian distribution with 

Pa{Ji,j) = exp(-j5/2J 2 )/(J^) (50) 

and (iii) the Laplacian distribution with 

P L (Jr,j) = exp(-\/2 | J itj /J \)/(jV2). (51) 

All distributions are symmetric with zero mean and variance (Jfj) / J 2 = L The 
dynamics of the model is given by a master equation where the rates are chosen 
according to heat-bath dynamics. It is by now established, see e.g. [82] for a review, 
that for d > 2 dimensions this model undergoes an equilibrium phase-transition 
between a paramagnetic and a frustrated spin-glass phase. There has been consider- 
able debate on the precise relationship between the relevant time and length scales 
for quenches below the spin-glass critical temperature T c . It has been attempted to 
summarize the present state of knowledge into the form [19] 

m ~ L '^{Hmf) (52) 

where Ao is an energy scale of order T c , tp is a barrier exponent and £(T) is the 
equilibrium correlation length at temperature T. This form has been used to fit 
successfully simulational data in the three-dimensional and in the four- dimensional 
Edwards- Anderson model [19, 18]. Although the typical length scales are merely of 
the order of a few lattice sizes, see e.g. [121], the relaxation times are sufficiently large 
for a dynamical scaling to set in. While the expression (52), if correct, points towards 
a cross-over behaviour between a simple power-law scaling and an exponential scaling 
for T < T c in spin glasses as would follow from the droplet model, it also suggests 
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Fig. 7. Scaling of (a) the autocorrelation function C(ys,s) and (b) the thermore- 
manent magnetization M(ys, s) of the critical three-dimensional Edwards- Anderson 
spin glass with a bimodal distribution of the couplings for (from top to bottom) 
y — 5, 8, 10, 15. The full curves are C(ys, s) = cos~ b and M(ys, s) = mos~ a , where 
Co, mo were fitted to the numerical data, (c) Temporal evolution of the autocorrela- 
tion C(t, 0) in the three-dimensional cases (b: bimodal, g: Gaussian, 1: Laplacian). 



that at criticality simple power-law scaling (in the sense discussed above in section 3) 
should prevail. 

From now on we consider the three- and the four-dimensional Ising spin glass 
quenched to T = T c from a fully disordered initial state and discuss the dynamical 
scaling behaviour of critical two-time quantities. We consider two classes of observ- 
ables: (a) thermoremanent magnetization where h is constant up to a single jump 
and (b) alternating susceptibility, where h = h(t) is oscillatory. 

4.1 Thermoremanent magnetization 

Let us first verify that the standard scaling forms (5) and (6) for the autocorrelation 
and the autoresponse are indeed observed in critical spin glasses. From the data 
shown in figure 7 for the case of a bimodal distribution of the couplings in three 
dimensions we conclude that both the correlation (figure 7a) and integrated response 
(figure 7b) are consistent with a power-law scaling. Similar results are obtained 
[63, 105] for the two other distributions, both in three and four space dimensions. 

From the slopes of figure 7 we obtain the estimates for the exponents a and 
b gathered in table 6. Interestingly, we observe that a — b within our numerical 
errors for a fixed choice of the distribution of the couplings, but that the results 
for different distributions are different. Even more intriguingly, other critical non- 
equilibrium dynamical quantities, as for example the exponent Xc/z describing the 
decay of the correlations in the long-time limit or the limit value of the fluctuation- 
dissipation ratio (8), also vary strongly and systematically with the form of the 
interaction distribution, see table 6 and figure 7c [105]. Here we emphasize the 
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Table 6. Nonequilibrium quantities of the critical Ising spin glass for bimodal, 
Gaussian, and Laplacian distributions of the nearest-neighbour couplings in d — 3 
and d = 4 dimensions. 





bimodal 


d 


a 6 Ac/z -Xoo 


3 
4 


0.060(4) 0.056(3) 0.362(5) 0.12(1) 
0.18(1) 0.180(5) 0.615(10) 0.20(1) 




Gaussian 


d 


a b Xc/z Xoo 


3 
4 


0.044(1) 0.043(1) 0.320(5) 0.09(1) 
0.169(4) 0.171(2) 0.58(1) 0.175(10) 




Laplacian 


d 


a b Xc/z Xoo 


3 
4 


0.033(3) 0.032(2) 0.259(2) 0.055(2) 
0.143(5) 0.140(3) 0.54(1) 0.13(1) 



diagnostic usefulness of the universal limit fluctuation-dissipation ratio Xoc whose 
numerical estimates vary up to a factor of 2 between the various distributions. 
Having excluded the most probable sources of systematic errors [105], we interpret 
the numerical data as strong evidence that in spin glasses critical out-of-equilibrium 
quantities do depend on the exact form of the distribution of the couplings. 9 It 
is not yet clear whether this intriguing property is only a purely dynamical effect. 
Whereas in the past some studies of static critical quantities observed a similar 
dependence on the distribution of the couplings [15], recent studies yielded for the 
different distributions values which agree within error bars [79, 81]. 

Let us now turn to the scaling functions themselves. Plotting the rescaled auto- 
correlator s C(t,s) versus the scaling variable t/s yields a nice data collapse com- 
patible with a simple power-law scaling L(t) ~ t 1 '" [63, 105]. In figure 8 we show 
the two-time scaling of the integrated response where again a collapse of the data 
in terms of a simple power-law scaling is observed. The fact that both the autocor- 
relator and the thermoremanent magnetization can be described in terms of such a 
power- law scaling is evidence in favour of the time-dependent length-scale (52). In 
principle, one would like to extract an exponent X' R /z from the slopes of the ther- 
moremanent magnetization. It turns out that the values of X' R / z thus obtained are 
significantly different from the ones found for Xc/z. For example, for the Ising spin 
glass with a bimodal distribution of the couplings one finds X' R /z = 0.45 (resp. 0.72) 
in three resp. four dimensions. On the other hand, since Xoo is finite, one must have 
Xr = Ac, but it may be necessary to go to very large values of y — t/s in order to 
see this. Indeed, if one considers larger values of y = t/s as is shown in the insets of 
figure 8ab, one observes a passage from an effective exponent X' R /z at intermediate 
values of the scaling variable y to the truly asymptotic value X R /z at larger values 
of y. Because of this passage from an effective exponent X' R /z at intermediate values 

9 We recall that a similar result was already seen for the disordered, but unfrus- 
trated, Ising model in section 3. 
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Fig. 8. Scaling of the thermoremanent magnetization in the Edwards- Anderson spin 
glass at criticality with a bimodal distribution of the couplings in (a) three and (b) 
four dimensions. The full curve is the prediction (43) of local scale-invariance. The 
results of very longs runs for a single waiting time are shown in the insets where (a) 
s — 100 in three dimenions and (b) s = 25 in four dimensions. 



of y = t/s to the truly asymptotic value for larger y we cannot expect the LSI 
equation (43) to describe the scaling function /m(j/) for all values of y = t/s. In 
figure 8ab we compare the numerical data, in both three and four dimensions with 
the prediction (43) where we have inserted the values of the exponents a and X'r/z. 
We find a nice agreement of the prediction (43) of local scale-invariance with our 
data for y = t/s not too large, but for very large arguments the precise behaviour 
of the scaling function /m(j/) cannot be fully reproduced. 



4.2 Alternating susceptibility 

An interesting possibility to further test the behaviour of critical spin-glass models 
consists of working with a time-dependent (oscillating) magnetic field and to study 
simultaneously the dependence on time and on the imposed oscillation angular fre- 
quency u. For a harmonic magnetic field, it is common to consider the real and the 
imaginary part of the magnetic susceptibility 

x'W) 

where R(t,s) is the linear response discussed above. In this setting, 1/uj provides 
the second time-scale, the natural scaling variable is y = uit and the scaling regime 
should be reached in the limit uj — > and t — » oo. In many experiments and simu- 
lations, one averages over at least one period of the oscillating field, see e.g. [101]. 



= / du R(t, u) cos (co(t — it)) 



du R(t, u) sin (oj(t — it)) 



(53) 
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Then in a great variety of glass-forming substances quenched to below or near to 
their glass transition point one observes good but not always perfect evidence for 
an u;i-scaling behaviour of the following form of the period-averaged dissipative 
(imaginary) part [78, 117] 

X 77 ^, t) = XstH + XagcO, *) , XagcO, t) ~ K ge {iOty b " (54) 

where x"t is thought of as a 'stationary' contribution while the ageing behaviour is 
described by x"gc- The amplitude A" gc and the exponent b" are obtained from fits to 
the experimental data. Similar scaling forms have been proposed for the dispersive 
(real) part x' but in practice the imaginary part is usually easier to measure. It is 
usually thought that b' = b" . 

Assuming the validity of LSI, we have derived the scaling relation [62] 

b' = b" = a — a (55) 

where a' is the exponent defined in eq. (34). Hence this method gives direct access 
to the exponent a — a' and should be useful for future direct tests of whether a and 
a' are really different. 

In order to derive (55), we notice the central role played by the time difference 
r = t — u in (54) [56]. Depending on its value either an equilibrium behaviour or 
else an ageing behaviour is obtained. Recall from section 3 that there is a time-scale 
t p ~ r with < £ < 1 on which the transition between the two regimes occurs [122] 
such that R(t,s) ~ R cq (t — s) for t — s < t p and R(t,s) = s _1_a /jj(t/ 's) as given 
in eq. (34) for t — s > t p . In addition, one measures for u w t the response with 
respect to a change in the initial conditions and then instead of (34) one expects 
R ~ Rini{t) ~ t~ XR / z [20]. 10 . We therefore must introduce a further time-scale t s 
such that t — t e — O(l). Splitting the integral into three terms corresponding to 
these three regimes, one has 

x"(w,t)= / dr R(t, t — t) sincjr 
Jo 

rt p rt e 
= / dr R(t,t — r) sin u;r + / dr R(t, t — r) sin lot 
Jo J t p 

+ / dr R(t, t — r) sincjr 

Jt e 

~ / dr R C q(r) sin lot + t~ a / dv f R I 
Jo Jtp/t V 

I 

-xfM+.-jr'd./.^^+o^.) 

= x"H + *"°Xa'(wt) + O (t- x " /z ) (56) 

In the third line, we used the asymptotic forms of R(t, s) as described above. This 
means that the cross-over between the equilibrium and the ageing regimes is assumed 



smutty 



Jtp/t 

f drcosinwr 



t E /t 

in f„ I . 

1 - v J (1 - v) 1+a 



A similar argument is used in the derivation of eqs. (42,43) [56, 60] 
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to be very rapid. While this is certainly correct for simple magnets, its validity in 
spin-glasses is far from obvious. In the last line, we restricted ourselves to the long- 
time limit t — > oo. Here, the function x"( w ) merely depends on the equilibrium 
form of the linear response R cq (t,s). It also becomes clear that the often-found 
stationary term in the integrated response [29, 31, 34, 40, 78, 82, 108, 109, 117] does 
not necessarily require the separation of a similar 'stationary' part in the response 
function R(t,s) itself. An analogous expression can be derived for x'(w,t). 

We now insert the prediction (34) of LSI. Then the scaling functions read, to- 
gether with their leading behaviour as y — > oo [62] 

x »(y) = f B(l-a',^-a} y 1 -* 

^ f"y a '- a + f2'y- XR/z sin (y + I [o - Xr/z] 
Xziv) = foB ^-o', ^j- - 

X2F H~ '^ _ '2' — 2- + 2l' 2 + 2l'"Tj (58) 

* /i/~ a + /22/" Ah/z cos (y+\[a- Xr/z}) 

where 2-F3 is a hypergeometric function and f[ 2 and /" 2 are known constants pro- 
portional to the normalization constant fo. We see that there appear terms which 
decrease monotonously with y but that there are also oscillating terms. They are 
described by different exponents and must be extracted by a different experimental 
setup. The oscillating terms follow the oscillations of the external field and the de- 
crease of the oscillation amplitude gives a direct access to the exponent Xr/z. On the 
other hand, in many experiments the data are averaged over one or several periods 
of the external field. For y sufficiently large, the contribution of the oscillating term 
in eqs. (57,58) vanishes after averaging and then only a simple algebraic component 
remains, which permits to extract the exponent a — a! . For period-averaged data or 
else if Xr/z > a — a' , the leading behaviour for large arguments is 

X2{y) ~ xUy) ~ y'- a (59) 

and the scaling relations (55) follow. Since the new exponent a! enters in (55), this 
suggests that no scaling relation between b' = b" and the other ageing exponents, 
with a' excluded, exists. 

We now describe a test [62] of the predictions eqs. (57,58) in the three- 
dimensional Ising spin glass with bimodal disorder, quenched to its critical point 
T c w 1.19. In order to study the alternating susceptibility far from equilibrium we 
prepared the system in an uncorrelated initial state (corresponding to an infinite 
initial temperature) before quenching it to the critical point at time t — 0. At the 
same time an external oscillating magnetic field 

h(t) — ho cos cot (60) 

was switched on, with its amplitude fixed at ho = 0.05 which is well inside the linear- 
response regime. We consider different values of the angular frequency u> = 2tv /p 
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with p ranging from 50 to 1600. Numerically, the in-phase and the out-of-phase 
susceptibilities are given by the expressions [4] 

x" (u,t) = m(t) sin cut , x ( w > t) = m(t) cos cut (61) 

with m(t) = o"i(t) but the equilibrium parts of %" an d x' must be subtracted off 

[62]. 
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Fig. 9. Scaling of the dissipative part y^iut) of the alternating susceptibility as 
function of the scaling variable uit for different angular frequencies uj — 2ir/p with 
p — 1600, 800, and 400. The full curve is the theoretical prediction (57) with fo = 
0.002 and a' = —0.70 but which has also been shifted horizontally by y — > y + Ay, 
with Ay — —0.45, see text. Statistical error bars are smaller than the symbol sizes. 




In figures 9 and 10 [62] we compare the expected scaling behaviour of the ageing 
part 

X2 = X2M) and x'2 = X2M). (62) 
For the larger values of p, corresponding to the smaller values of u>, one observes a 
very good data collapse for both quantities, which furnishes clear evidence in favour 
of a power- law scaling at T = T c . The value a — 0.060(4) determined previously 
from the decay of the thermoremanent magnetization [63, 105] was used. For smaller 
values of p, the collapse is less good. It is possible that for the corresponding values 
of uj the dynamical scaling regime is not yet reached. 

The data shown in these two figures can in principle be compared directly with 
the analytical predictions (57) and (58), provided however, that the positions of 
maxima of x'2 an< i X2 are shifted by the amount Ay w —0.45 when compared with 
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Fig. 10. The same as in Figure 9, but now for the dispersive part )(^{<jJt). The full 
curve is the shifted theoretical prediction (58) with /o = 0.002 and a' = —0.70. 



the positions obtained in the analytical treatment of the previous section. Where 
this shift comes from is not understood. It is possible, however, that the cross- 
over between the equilibrium and the ageing regimes is not almost instantaneous in 
contrast to what we assumed in the derivation of eq. (56). 

Since both a and Xr/z were known before [63, 105] the only free parameters in 
this comparison are the amplitude /o and the exponent a' . The final estimates are 
[62] 

a = -0.70(3) , f = 0.00203(1) (63) 

and the fit is compared to the data in the figures. These values (63) of the parameters 
describe consistently both x'2 and X2- Although some discrepancies are observed for 
small values of y = uit, the overall agreement between the simulation and the shifted 
theoretical prediction is very good. Clearly, LSI captures an essential part of the 
ageing behaviour of critical spin glass models. 

We finish with a discussion of existing experimental results on the scaling of the 
alternating susceptibility. A detailed discussion of a possible scaling of the alternat- 
ing susceptibility was presented by Suzuki and Suzuki for the short-ranged Ising spin 
glass Cuo.sCoo.sC^-FeChj graphite bi-intercalation compound [117]. After a rapid 
quench to below the glass temperature T g = 3.92(11)K, they measure x"( w ii) f° r 
fixed lo and find that their (period-averaged) data are well-fitted by the power-law 

FKt) = XoH+^>)*- 6 " (64) 

While the fitted exponent b" depends only slightly on w, they further show evidence 
for a power-law A" (u) = A'quj~^ and find that . . the value of //' is almost the 
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same as that ofb"" [117]. In this way, they arrive at the wt-scaling form 

F(w,*) = xi'(w)+^'(w*)" fc " (65) 

and a similar form for x'( w j t) where ". . . b' and b" are of the same order at the same 
temperature" [117]. Experimentally measured values of the exponent b" (and also 
b') of some materials are collected in table 7. Comparing the experimental results 
of Suzuki and Suzuki [117] with the theoretical scaling form (56), one sees: 

1. the experimental evidence for a pure cjt-scaling indicates that the exponent a 
must indeed be very small. 

2. when considering the leading behaviour for y = cot large (their data go up to 
V IS 10 6 [H7]) and recalling that the experimental data are averaged over at 
least one period of the external field, one can read off 

b" — a — a , n" — —a (66) 

and the observed [117] near equality b" ~ /i" is again consistent with a being 
negligibly small. 

3. the data are consistent with the theoretically requested relation b' = b" . 



Table 7. Measured values of the exponents b" and b' in several glassy materials, us- 
ing the scaling form (54). Here T g stands for the glass transition temperature and T is 
the temperature where the data were taken. For Feo.sMno.sTiOs and CdCr1.7Ino.3S4 
the relation b' — b" was assumed. The simulational results in the critical Ising spin 
glass are also included. 



Material 


T g [K] T[K] b" b' 


Ref. 


Cuo.sCoo.sCb-FeCla 
- GBIC 


3.92(11) 3.25 0.01(4) 0.08(3) 
3.5 0.017(32) 0.05(2) 
3.75 0.16(3) 0.20(2) 
3.85 0.15(3) 
3.95 0.16(4) 0.20(2) 


[117] Ising 

spin glass 


Feo.sMno.sTiOa 


20.7 15 0.14(3) 
19 0.14(3) 


[39, 37] Ising 

spin glass 


CdCr1.7Ino.3S4 


16.7 12 0.18(3) 
14 0.18(3) 


[39, 37] Heisenberg 
spin glass 


CdCr 2l In2-2^S4 x = 0.95 
x = 0.90 


70 8 0.2 
67 0.2 
50 42 0.20 


[38, 37] disordered 
ferromagnet 

[119, 37] 


Pb(Mgi /3 Nb2 / 3)03 


~ 220 < 220 0.17 


[31] relaxor 

ferroelectric 


theory (T = T c ) 


0.76(3) 0.76(3) 


[62] Ising 

spin glass 



Similar values of b" were observed for several other materials, quite indepen- 
dently of the precise physical nature as can be seen from table 7, but the errors are 
still too large to permit a discussion of the universality of the exponents. However, 
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the experimental data are in many of these systems at least as well described by a 
logarithmic scaling as expected from the droplet theory [37, 78, 39]. Furthermore, 
in several systems also strong deviations from a simple c<;t-scaling were found, see 
[32]. Finally, we mention that in systems like /3-hydroquinol-clathrate [108] or even 
simple liquids like glycerol [85] a power-law dependence of the form xl'gc ~ t~ a or 
Xago ~ t~ a was observed. All in all, it is at present not completely clear why one 
should find a simple at-scaling of the alternating susceptibility whereas subageing 
is frequenlt admitted for the thermoremanent magnetization. 

A last important observation follows from comparing the experimental data for 
the exponents b' = b" ~ 0.1 — 0.2 for all the 3D materials studied so far with the 
theoretical estimate b" ~ 0.76 derived from the 3D Ising spin glass with binary 
disorder. The values are very far from each other and this huge difference calls 
for an explanation. Could this be seen as an indication that the spin glass models 
considered by theorists only capture imperfectly what is going on in real materials ? 



5 Discussion 

The available evidence for universality in the ageing behaviour seen in either glassy 
materials, disordered systems or else simple many-body systems without disorder 
may be taken, in combination with dynamical scaling, as suggestive of the existence 
of deeper dynamical symmetries in such systems. We have advocated here the point 
of view that such a new symmetry might be sought through an extension of standard 
dynamical scaling towards a local scale-invariance. The present formulation of the 
theory has been built around dynamical scaling as it occurs in full ageing, but 
this was mainly for reasons of technical simplicity and one could consider different 
forms of dynamical scaling, with a modified form of the Lie algebra generator Xo, 
if required. 

Even within this specific context, it is in principle possible to generalize in many 
different ways. For that reason, it is important to be able to test quantitatively some 
of the consequences of local scale-invariance. Here the calculation of the two-time 
response function R(t, s) is one of the easiest tasks and has the advantage that 
prediction are relatively easy to test in simulations, for example by considering an 
integrated response such as the thermoremanent magnetization Mtrm((, s). 

The central idea of LSI is the decomposition of the stochastic Langevin equation 
into a deterministic part which can possess local scale-invariance and into a noise 
part which does break this. Remarkably, the noise part is usually built in such a 
way that any n-point function of the full theory can be reduced exactly to a certain 
m-point function of the noiseless theory (where m > n). This reduction depends 
on the Bargman superselection rules which are proven only for z — 2 at the time 
of writing. Then the symmetries of the noiseless part yield conditions on the form 
of the m-point functions. In certain cases, such as the two-time response functions, 
these constraints are enough to fix the functional form, see eq. (34). 

In models where the deterministic part of the Langevin equation is linear, LSI 
can be proven, for example in the spherical model [103], the ID Glauber-Ising model 
[103, 65], bosonic particle-reaction models [9] or certain simple growth models [110, 
12]. Furthermore, there is a systematic way to construct LSI- invariant non- linear 
equations [115, 9]. Since these are not identical to the Langevin equations one usually 
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starts from [71], numerical evidence for those systems is crucial, in the absence of 
exact solutions. 

In order to be able to carry out meaningful tests, one must first show that simple 
ageing (on which LSI is built) holds true. For non-conserved ferromagnets this is 
thought to be certain [20] and we have discussed the evidence that this remains so 
for disordered, non-glassy Ising models. In most numerically studied systems, only 
Mtrm(£, s) has been compared with LSI. Provided that corrections to scaling are 
treated carefully, in all cases the results are consistent with LSI, see tables 2 and 
3. The fact that Mtkm (£, s) as measured in a noisy system is consistent with the 
integral of (34) found from the symmetries of the noiseless theory, even if z 2, is 
a strong indication that an extension of the Bargman superselection rules to z 7^ 2 
should exist. 

At present, the only other observables measured have been space-time responses 
and autocorrelations in the two- and three-dimensional Ising models[58, 61] and 
in the two-dimensional Potts-g models (q = 2,3,8) [89]. It would be important to 
consider different quantities, e.g. alternating susceptibilities. 

When turning to critical spin glasses, although there is evidence that simple 
ageing applies, the agreement of the data with LSI was only partial. It should be 
recalled that the analysis performed relied on a rapid change between quasi-static 
and scaling behaviour, which might be an oversimplification. Remarkably, in the 
presence of disorder we have seen evidence that the non-equilibrium exponents may 
depend on the distribution of the random exchange couplings Jjj, independently 
on whether the system was glassy or not. Since this is a highly controversial issue, 
further tests would be most welcome. 

At present, the available evidence 11 from both analytical solved models and nu- 
merical studies supports the predictions of local scale-invariance for the two-time 
response (and in some cases also correlation) functions. This suggests that the kind 
of Langevin equation used in the description of these models should possess a funda- 
mental, hitherto unsuspected and non-trivial dynamical symmetry underlying cer- 
tain stochastic Langevin equations. Already the present examples show that the cal- 
culation of time-dependent quantities such as the thermoremanent magnetization is 
considerably simplified and leads to results not available as yet by other methods. 
Local scale-invariance, similar in spirit to the tremendously successful techniques 
of conformal field-theory, has the potential of leading to a much more profound 
understanding of non-equilibrium critical phenomena. 

Finally, comparison of simulational results with experimental evidence for both 
glassy and non-glassy systems (see tables 5 & 7) points to some systematic differ- 
ences which need to be clarified. 
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11 Of course models such as the celebrated zero-range process, see [42, 50], which 
lacks a proper spatial structure should be excluded from this kind of consideration. 
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